Level Crossing Analysis of Growing surfaces 
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We investigate the average frequency of positive slope u„, crossing the height a — h — h in the 
surface growing processes. The exact level crossing analysis of the random deposition model and the 
Kardar-Parisi- Zhang equation in the strong coupling limit before creation of singularities are given. 
PACS: 52.75.Rx, 68.35.Ct. 



I. INTRODUCTION 

Due to the technical importance and fundamental in- 
terest, a great deal of effort has been devoted to under- 
standing the mechanism of thin-film growth and the ki- 
netic roughening of growing surfaces in various growth 
techniques. Analytical and numerical treatments of sim- 
ple growth models suggest, quite generally, the height 
fluctuations have a self-similar character and their aver- 
age correlations exhibit a dynamic scaling form [1-6]. It 
is known that to derive the quantitative information of 
the surface morphology one may consider a sample of size 
L and define the mean height of growing film h and its 
roughness w by [1]: 



L/2 



L/2 



dxh(x, t) 



and 



w(L,t) = (((h-h) 2 )) 



1/2 



(1) 



(2) 



where (• • •) denotes an averaging over different realiza- 
tion (samples) . Starting from a flat interface (one of the 
possible initial conditions), it was conjectured by Family 
and Vicsek [7] that a scaling of space by factor b and 
of time by a factor b z [z is the dynamical scaling expo- 
nent), re-scales the roughness w by factor b x as follows: 
w(bL,b z t) = b x w(L,t), which implies that 



w{L,t)=L*f{t/L z ). 



(3) 



If for large t and fixed L (t/L z — > oo), w saturate then 
f(x) — > const., as x — ► oo. However, for fixed large L 
and 1 << t << L z , one expects that correlations of the 
height fluctuations are set up only within a distance t x / z 
and thus must be independent of L. This implies that for 
x << 1, f(x) ~ x 13 with [3 = x/ z. Thus dynamic scaling 
postulates that, w(L,t) ~ t 13 for 1 << t << L z and 
~ L x for t >> L z . The roughness exponent x an d the 



dynamic exponent z characterize the sclf-affinc geometry 
of the surface and its dynamics, respectively. 

Here we introduce the level crossing analysis in the 
context of surface growth processes. In the level cross- 
ing analysis we are interested in determining the average 
frequency ( in spatial dimension ) of observing of the def- 
inite value for height function h — h = a in growing thin 
films, from which one can find the averaged number 
of crossing the given height in sample with size L. The 
average number of visiting the height h — h = a with pos- 
itive slop will be N+ = v+L. It can be shown that the 
v+ can be written in terms of joint probability distribu- 
tion function (PDF) of h — h and its gradient. Therefore 
the quantity v+ carry the whole information of surface 
which lies in joint PDF of height and its gradient fluctu- 
ations. This work aims to study the frequency of positive 
slope crossing (i.e. in time t on the growing surface 
in a sample with size L. We introduce a quantity N^~ ot 
which is defined as N^ ot = v^da to measure the 

total number of crossing the surface with positive slop. 
The Nfo t and the path which is constructed by growing 
surface are in the same order. It is expected that in the 
stationary state the N^ ot to become size dependent. We 
determine the time and height dependence of v£ for two 
exactly solvable models, random deposition model (RD) 
and Kardar-Parisi-Zhang (KPZ) equation in the strong 
coupling limit and before creation of singularities (sharp 
valleys) with short range forcing. It is shown that the 
RD-model and KPZ equation have different i/+ but N^ ot 
scales with time as t x l 2 in both models. 

In section 2 we discuss the connection between v+ and 
underlying probability distribution functions (PDF) of 
growing surfaces. Exact expression of j/+ for the RD 
model and with short-range forcing is given in section 
3. In section 4 we derive the integral representation of 
for the KPZ equation in 1+1 dimension and in the 
strong coupling limit before the creation of singularities. 
We summarize the results in section 5. 



II. THE LEVEL CROSSING ANALYSIS OF 
GROWING SURFACE 



h(4 



Consider a sample function of an ensemble of functions 
which make up the homogeneous random process h{x,t). 
Let n+ denote the number of positive slope crossing of 
h{x) — h = a in time t for a typical sample size L (see 
figure 1 ) and let the mean value for all the samples be 
N+(L) where 



N+(L)=E[n+(L)}. 



(4) 



Since the process is homogeneous, if we take a second 
interval of L immediately following the first we shall ob- 
tain the same result, and for the two intervals together 
we shall therefore obtain 



N+(2L) = 2N+(L), 



(5) 



from which it follows that, for a homogeneous process, 
the average number of crossing is proportional to the 
space interval L. Hence 



or 



N+(L) cx L, 



N+{L) = V +L. 



(6) 



(7) 



which is the average frequency of positive slope cross- 
ing of the level h — h = a. We now consider how the 
frequency parameter can be deduced from the under- 
lying probability distributions for h— h. Consider a small 
length dl of a typical sample function. Since we are as- 
suming that the process h — h is a smooth function of 
x, with no sudden ups and downs, if dl is small enough, 
the sample can only cross h—h — a with positive slope 
if h — h < a at the beginning of the interval location x. 
Furthermore there is a minimum slope at position x if the 
level h — h = a is to be crossed in interval dl depending 
on the value of h — h at location x. So there will be a 
positive crossing of h — h = a in the next space interval 
dl if, at position x, 



h— h < a 



and 



d(h - h) > a-(h-h) ^ 



dl 



dl 



Actually what we really mean is that there will be high 
probability of a crossing in interval dl if these conditions 
are satisfied [8,9]. 

In order to determine whether the above conditions are 
satisfied at any arbitrary location x, we must find how 
the values of y — h — h and y' = ^ are distributed by 
considering their joint probability density p(y,y'). Sup- 
pose that the level y = a and interval dl are specified. 
Then we are only interested in values of y < a and values 
of y' = (^) > , which means that the region be- 
tween the lines y = a and y' = ^-rf^ in the plane (y.y'). 





h-h=a 



FIG. 1. positive slope crossing of the level h — h — a. 

Hence the probability of positive slope crossing of y = a 
in dl is: 



/ dy' / dyp(y,y'). 

JO Ja — y'dl 



When c/Z — > 0, it is legitimate to put 

p(y,y') =p(y = a,y') 



(9) 



(10) 



Since at large values of y and y' the probability density 
function approaches zero fast enough, therefore eq.(6) 
may be written as: 



dy' 



dyp(y = a, y') 



(11) 



a—y'dl 



in which the integrand is no longer a function of y so 
that the first integral is just: j^-y'di^yPiy = a >y') = 
p(y = a,y')y'dl, so that the probability of slope crossing 
of y — a in dl is equal to: 



dl / p{a,y')y'dy' 
Jo 



(12) 



in which the term p(a, y') is the joint probability density 
p{y,y') evaluated at y = a. 

We have said that the average number of positive slope 
crossing in scale L is v£L, according to (7). The average 
number of crossing in interval dl if therefore v+dl. So 
average number of positive crossings of y = a in interval 
dl is equal to the probability of positive crossing of y = 
a in dl, which is only true because dl is small and the 
process y(x) is smooth so that there cannot be more than 
one crossing of y = a in space interval dl, Therefore we 
have v^dl = dl p(a,y')y'dy' , from which we get the 
following result for the frequency parameter in terms 
of the joint probability density function p(y, y') 



f 

Jo 



p{a,y')y'c 



(13) 



In the following sections we are going to derive the v£ 
via the joint PDF of h— h and height gradient. To derive 
the joint PDF we use the master equation method [10-11]. 
This method enables us find the v+ in terms of generating 
function Z(X, u, x, t) — (exp(—i\(h(x,t) — h)—ifiu(x,t))), 
where u(x,t) = — Vh. 



III. THE FREQUENCY OF A DEFINITE 
HEIGHT WITH POSITIVE SLOPE FOR THE 
RANDOM DEPOSITION MODEL 

In random deposition model particles are dropped ran- 
domly over deposition sites, and stick to the top of the 
pre-existing column on the site [ 1 ]. The height of each 
column thus performs an independent random walk. This 
model leads to unrealistically rough surface whose over- 
all width growth with the exponent (3 = \ without sat- 
uration. In the continuum limit the random deposition 
model is described by the following equations; 



d_ 

dt 



h(x,t) = f(x,t) 



d_ 

dt 



u(x, t) = f x 



(14) 



where h(x,t) is the height field, u{x,t) = J^h(x,t) and 
f(x,t) is a zero mean random force gaussian correlated 
in space and white in time, 



(f(x, t)f(x', t')) = 2D D(x - x')6(t - t') 



(15) 



where D(x) is space correlation function and is an even 
function of its argument. It has the following form; 



D(x - x') = 



y/ircr 



exp(-^-#) (16) 



where a is the standard deviation of D(x — x'). The 
average force on the interface is not essential and can 
be removed by a simple shift of h ( i.e h — ► (h — Ft) 
where F =< f(x,t) >) . Typically the correlation of 
forcing is considered as delta function for mimicking the 
short range correlation. We regularize the delta function 
correlation by a gaussian function. When the variance a 
is much less than the system size we would expect that 
the model would represent a short range character for 
the forcing. So we would stress that our calculations are 
done for finite a <C L, where L is the system size. The 
parameters Dq is describing the noise strength. 

Now assuming the homogeneity we define the generat- 
ing function as: 



Z(X 7 fi, t) = (cxp(—iXh(x, t) — i[iu(x, t)) 



(17) 



by using the eq.(14) we can find the following equation 
for the evolution of Z(X, fi, t): 



d_ 

dt 



Z(X, \i, t) = —iX(h t (x, t) exp(—iXh(x, t) — i[iu(x, t)) 
—i/j,(u t (x, t) exp(—iXh(x, i) — i^iu(x, t)} 

= -i\(f(x, t) cxp(-i\h(x, t) - i[ij x {x, t)) 
—i/j,(ut exp(—i\h(x, t) — ifiu(x, t)) 

= -X 2 D D(0)Z + ii 2 D D xx (0)Z. 



The joint probability density function of h and u can 
be obtain by fourier transform of the generating function: 

P{h,u,t) = -!- I ' dXdne iXh+iilu Z{X,n,t), (19) 

Z7T J 

so by fourier transforming of the eq.(18) we get the 
Fokkcr-Planck equation as: 



d 2 

D D xx (0)—P. (20) 



the solution of the above equation can be separated as 
P(h,u,t) = pi(h,t)p2(u,t). Using the initial conditions 
Pi(/i,0) = 5(h) and P 2 (u,0) = 5(u) (starting from flat 
surface) it can be shown that: 



P(h,u,t) 



1 



4irtV-D 2 D(0)D xx (0) 



exp(- 



h 2 



+ 



AD Q D(0)t AD Q D xx (0)t 



), (21) 



from which the frequency of repeating a definite height 
(h(x,t) — a) can be calculated as 



uP(a, u)du 



D xx (0) 



D(0) 



exp(- 



4L>(0)r 



4D Q D(0)t' 



(22) 



The quantity f+ in RD model has a gaussian form 
with respect to a. The zero level crossing scales with a 
as Va=a ~ cr _1 - Also using the eq.(22) it is found that 
N+ t = D^Tr-^a-VH 1 / 2 . This shows that there is no 
stationary state for the RD model and the quantity N^ ot 
diverges without saturation. 



IV. FREQUENCY OF A DEFINITE HEIGHT 
WITH POSITIVE SLOPE FOR KPZ EQUATION 
BEFORE THE SINGULARITY FORMATION 

In the Kardar-Parisi-Zhang (KPZ) model (e.g. in the 
1-dimension), the surface height h(x,t) on the top of lo- 
cation x of 1-dimensional substrate satisfies a stochastic 
random equation: 



^-^ dxh f=ud 2 x h + f(x,t). 



(23) 



where a > and / is a zero-mean, statistically homo- 
geneous, white in time and gaussian process with covari- 
ance as eq.(16). The parameters v, a and D a ( and a) 
are describing surface relaxation, lateral growth and the 
noise strength, respectively. Let us define the generating 
(18) function Z(X,u,x,t) as: 



Z(X, /j, x, t) = (exp(— i\(h(x, t) — h) — i/j,u(x, t))). 



Where u(x,t) — —d x h(x,t). Assuming statistical ho- 
mogeneity i.e. Z x = it follows from eq.(23) that Z 
satisfies the following equation; 

- i[iZ t = ^{t)\fj,Z - —XfiZ^ + iX 2 fik(0)Z 



- ii/k xx (0)Z - i{vX 2 + ia\)Z ll 

— /j, 2 v(u xx (x, t) exp(— iXh(x, t) — ifiu(x, t))). (24) 
where k(x-x') = 2D D(x-x'), j(t) = h t , k(0) = D ° 



and k xx (0) 



2-Do 



and h(x,t) — h(x,t) — h. 



Once 



trying to develop the statistical theory of the rough- 
ened surface it becomes clear that the inter-dependency 
of the height difference and height gradient statistics 
would be taken into account. The very existence of the 
non-linear term in the KPZ equation leads to develop- 
ment of the cusp singularities in a finite time and in the 
strong coupling limit i.e. v — > 0. So one would distin- 
guish between different time regimes. Recently it has 
been shown that starting from the flat interface the KPZ 
equation will develop sharp valleys singularity after time 
scale t*, where t* depends on the forcing properties as 
** = (i) 4/ W /6 AT 1/3 a~ 2/3 ff 5/3 [10]. This means that 
for time scales before t* the relaxation contribution tends 
to zero when v — > 0. In this regime one can observe that 
the generating function equation is closed. The solutions 
of the resulting equation is easily derived, ( starting from 
a flat surface, i.e. h(x,0) = and u(x,0) = 0), and has 
the following form [10]; 



Z(p, A, t) = (1 - tanh 2 ( v /2ifc xx (0)aAt)) 
exp[-- ln(l - t&xi\i A {yj2ik xx {0)a\t)) 



■ t&n\r 1 (ta,nh 2 {^2ik xx {Qi)aXt)) - X 2 k(0)t 




[...].. ..J. 



FIG. 2. plot of Va vs a for the KPZ equation in the strong 
coupling and befor the creation of sharp vallyes for time scale 
t/t* = 0.05,0.15 and 0.25. 



that the i/+ and N^~ ot can be written in terms of the gen- 
erating function Z as: 



2tt 



+00 /* + oo 



oo J — oo 



;Z(X, n) exp(iAa)<iA(i/i (26) 



1_ m 2 ( l - tanh(^2ik xx {0)a\t) 
16 n ( l + tanh(y/2ik xx (0)aXty 



and 



2 V a\ 

One can construct P(h, u, t) in terms of generating func- 
tion Z as eq.(19), from which the frequency of repeating 
a definite height (h(x,t) — h = a) with positive slope can 
be calculated as ^+ = J °° uP(a,u)du. 

In fig. 2 we plot the for time scales before creation 
of singularity, t/t* = 0.05,0.15 and 0.25. In the KPZ 
equation due to the nonlinear term there is no h — > —h 
symmetry and one can deduce that the ^+ is also not 
symmetric under h — > —h. 

To derive the N^~ ot let us express the v+ and N^ ot in 
terms of the generating function Z. It can be easily shown 



, + oo 1 

/ -Z(X^Q,^dX. 



(27) 



Using the eq.( 25 ) one finds N+ t ~ a^/H 1 / 2 . Wc 
note that the expression of the z/+ for the RD model and 
KPZ equation before t* are different functions of a but 
N^ ot scales as <~ er -3 / 2 ^/ 2 in both models. In fig. (3) 
using the direct numerical integration of joint PDF of 
height and its gradient we plot the N to t vs t. In this 
gragh Ntot ~ t 1 / 2 which is in agreement with analytical 
prediction. 




FIG. 3. log-log plot of Ntot (A) vs t for the KPZ equation 
in the strong coupling before the creation of sharp vallyes. 
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V. CONCLUSION 



We obtained some results in the problems of RD model 
and KPZ equation in 1+1 dimensions with a Gaussian 
forcing which is white in time and short range correlated 
in space. We determined the explicit expression of av- 
erage frequency of crossing i.e. z/ + of observing of the 
definite value for height function h — h = a in a growing 
thin films for the RD model, from which one can find the 
averaged number of crossing the given height in a sam- 
ple with size L. It is shown the v+ is symmetric under 
h — > —h. The integral representation of ^+ is given for 
the KPZ equation in the strong coupling limit before the 
creation of sharp valleys. We introduced the quantity 
Ntot — f-tt v a^ a -> which measures the total number of 
positive crossing of growing surface and show that for the 
RD-model and the KPZ equation in the strong coupling 
limit and before the creation of sharp valleys N^ ot scales 
as cr -3 / 2 * 1 / 2 . It is noted that for these models v+ has 
different expression in terms of a. The ideas presented 
in this paper can be used to find the v + of the general 
Langevin equation with given drift and diffusion coeffi- 
cients. 
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